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ABSTRACT 

We have constructed a photoionization model, using the 3-D Monte Carlo code 
MOCASSIN for one of the hydrogen-deficient knots (J3) of the born-again planetary nebula 
C*") • Abell 30. The model consists of spherical knots, comprising a cold, dense, hydrogen-deficient 

core with very high metal abundances. The inner core, occupying 9.1% of the total volume 
of the knot, is surrounded by a less dense hydrogen-deficient and metal-enriched gas enve- 
lope, with less extreme abundances. The envelope of the knot might have been formed by the 
^-^ \ mixing of the knot material with the surrounding nebular gas. 

This bi-chemistry, bi-density model did not produce enough heating to match the fluxes of 

■ the collisionally excited emission lines (CELs) and of the optical recombination lines (ORLs) 

observed in the spectrum of the knot. We therefore included heating by photoelectric emission 

q ■ from dust grains in the thermal equilibrium calculations, and found that dust-to-gas ratios of 

?H ' 0.077 and 0.107 by mass for the central core and the envelope of the knot, respectively, are 

r/3 , sufficient to fit the spectrum. Surprisingly, photoelectric emission from grains is the dominant 

source of heating in the hot envelope of the knot, while heating by photoionization of helium 
and heavy elements dominates in the cold core. 

We obtain a good fit with the observations for most of the significant emission lines 
1 treated in our model. The two major discrepancies occurred for the [O II] 3727,29 A doublet 

J-j | and the [N II] 6548, 6584 A lines, which are severely underestimated in our model. Recombi- 

nation contributions could be significant and we included them for the O II transitions. How- 
ever, this was not sufficient to resolve the discrepancy, due to the high collisional de-excitation 
rates in the dense core, where most of the recombination lines would be produced. This possi- 
bly highlights a weakness using a discontinuous density distribution like ours, where in reality 
one might expect an intermediate phase to exist. 

The chemical abundances inferred from our modelling of the central core region and of 
the envelope of the knot ar e, at least qualitative ly, in agreement with the abundances derived 
by the empirical analysis of Wesso n et al J (12003). although the discrepancies between the core 
and the envelope abundances that we find are less dramatic than those implied by the ORL and 
CEL empirical analysis. Our models also indicate, in agreement with the empirical analysis 
of Wesson et al. (2003), that the C/O ratio in the two regions of t he knot is less t han unity, 
contrary to theoretical predictions for born-again nebulae (e.g. llben & Renzinill983l) . 
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1 INTRODUCTION 

The large angular diameter planetary nebula Abell 30 belongs to 
the very peculiar class of planetary nebulae with hydrogen-rich en- 
velopes surrounding hydrogen-deficient inner knots tTacob^ 1979^ 
- other examples include Abell 78 <Jacobvlll979h and Abell 58 
JPollacco etalJl992h . The outer envelope of Abell 30 has a radius 



of 64" and an almost perfectly round shape, showing a high excita- 
tion and a very low surface brightness. Given the size of the outer 
envelope and the fact that it appears to be very diffuse, Abell 30 is 
believe d to be quite old (^10,000 yr, i f a distance of 1.3 kpc is as- 
sumed; Harrington & Feibelman 1984). Its hydrogen-deficient cen- 
tral star, however, still appears to be quite bright (V=14.3). These 
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characteristics, and, in particular, the presence of fo ur symmetri- 
cally placed hydrogen-deficient inner knots led llben et alj 1 19831) 
to the classification of this object as a planetary nebula with a 
born-ag ain nucleus. Th e central knots w ere ide ntified simultane- 
ously bv ljacobvl 1 19791) and lHazard et alJ ^980); according to the 
born-again scenario, these knots consist of material whose abun- 
dance distribution has been altered by complete hydrogen burning 
and later ejected when the star returned to a brief red-giant phase 
after a final He-shell flash. The mechanism by which these knots 
would have been emitted is not known. The strongly axially sym- 
metric structure of Abell 30 led iLi viol ll995l) to suggest that tidal 
interaction of the central star with a companion may be a plausible 
mechanism. Following Jacoby's discovery, the knots were named 
Jl to J4, clockwise, starting from the southernmost one. 

Since the discovery of the hydrogen-deficient knots, Abell 30 
has been widely studied. Ijacobv & Ford ll983h obtained op- 
tical spectra and derived ab undances for knots J3 and J4. 
lHarrington & Feibelmanl ^984) obtained IUE spectra of the nebu- 
losity near the nucleus, in proximity to one of the knots, and carried 
out a study of the chemical abundanc es from the obse rved emis- 
sion lines. More recently. iGuerrero & Manc hado 1 19961) published 
a study of the chemical abundances of hydrogen-poor planetary 
nebulae, including Abell 30, based on long-slit data obtained in 
1986 at the 3 . 6 m te lescope at ESO, La Silla Observatory (Chile). 
IWesson et al1l2003l) have analyzed new observational data, includ- 
ing optical spectra of the polar knots J 1 and J3 and of the central star 
of Abell 30, obtained in 2000 with the ISIS spectrograph mounted 
on the 4.2 m WHT at the Observatorio del Roque de los Mucha- 
chos, on La Palma, Spain, as well as archive ultraviolet spectra 
of knots J3 and J4 obtained in 1994 with the Faint Object Spec- 
trograph aboard the Hubble Space Telescope. They found large 
discrepancies between the ionic abundances measured from opti- 
cal recombination lines (ORLs) and those measured from colli- 
sional ly excited lines (CELs). This had already been noticed by 
Ijacobv & Ford J 19831) who found a very high carbon abundance 
(almost half that of hydrogen) from ORLs, far higher t han that de- 
rived from CELs. Subsequently, Harrington & Feib elmar] 1 19841) 
proposed that this discrepancy could be explained if very large 
temperature variations were present within the knots of Abell 30. 
A high concentration of heavy elements in the core of the knots 
could create a very cool (~1000 K) but still ionized region emit- 
ting the observed ORLs. This could be surrounded by a hotter re- 
gion consisting of material with more typical nebular abundances 
emitting the observed CELs. The analysis of IWesson et alj 1200 3l) 
is in agreement with this interpretation. 

In Section 2 of this paper we give a description of the pho- 
toionization modelling, together with details of the parameters 
used. A strategy which was used in order to estimate the effects 
of photoelectric heating from dust grains is also outlined in this 
section. Our model results are presented in Section 3, while a dis- 
cussion of these and our final conclusions are given in Section 4. 



2 PHOTOIONIZATION MODELLING 
2.1 Preliminary modelling 

We have constructed a photoionization mode l for kno t J3 o f 
Abell 30 based on the observational data of IWesson et alj i2003l) . 
in order to investigate the temperature and ionic structures which 
could give rise to the observed ORLs and CELs. Given the very low 
surface brightness of the main extended nebula, and hence the very 



low density of material there, the effects of contamination of the 
spectra of the knots from this region are expected to be very small. 
Nevertheless, a background nebular spectrum wa s subtracted from 
the spectra of the knots by Wes son et alj |2003) in order to cor- 
rect for any contamination. The three-dimensional photoionization 
code, MOCASSIN, described bv lErcolano et alJ 120031) . was used in 
this work, since this code can self-consistently treat asymmetries 
and density and chemical inhomogeneities. Our model consisted 
of a three-dimensional Cartesian grid with the ionizing radiation 
source placed in one corner, with variable spatial sampling, with 
most of the cells being in the knot region, while the rest of the neb- 
ular region was assumed to be empty. 

The general procedure used to constrain the photoionization 
models descibed in this paper consists of an iterative process, in- 
volving the comparison of the predicted emission line fluxes, for 
some significant lines in the spectrum, with the values measured 
from the observations. The free parameters used in our models in- 
cluded nebular parameters, such as gas density and abundances, 
stellar parameters, such as luminosity, effective temperature and 
gravity, and, for the models described in Section l2~2l dust param- 
eters, such as dust-to-gas-ratio and grain size. The first step of the 
modelling process is to get a handle on the total gas mass and stel- 
lar luminosity required to match the total flux in the reference emis- 
sion line as well as the relative intensities of the strongest lines in 
the spectrum. In hydrogen-rich environments the H/3 line is gen- 
erally chosen as the reference. However, in the case of knot J3 of 
Abell 30 we have chosen the He I 5876 A line as the reference, 
since helium is the most abundant element. The exact match will, 
of course, also depend upon the ionization structure of the neb- 
ula, which is dependent on the gas density and also on the central 
star temperature. In general, all model parameters are dependent 
on each other and it is not possible to completely isolate the ef- 
fect of any given parameter on the resulting nebular structure and, 
consequently, on the nebular emission line spectrum. By looking 
at the diagnostic line ratios one can verify whether the temperature 
stucture of a given model is a realistic representation of the object 
being studied. The general ionization structure of the nebular gas 
is heavily dependent on the ionizing star parameters, as well as on 
the gas density distribution. The effective temperature and grav- 
ity of the central star were varied until the He 2+ /He + model ratio 
matched the value derived by the empirical analysis. The starting 
values were chosen to be L« = 5000 L^, T e ff = 1 10,000 K and 
log g = 6.0 lHarrington & Feibelman 1984). Hydrogen-rich non- 
LTE model atmospheres were used for our preliminary modelling, 
however for the later and final models we used the non-LTE hydro- 
gen deficient atmospheres described in Section l2~4l 

In the first model constructed, the density and abundances in 
the knot were assumed to be homogeneous. A first attempt was 
made by using a density of 3200 cm" 3 , a s inferred from the [O II] 
3726 A/3729A line ratio IWesson et ail EoQ3). Models were run 
both using the abundances derived from their empirical CEL analy- 
sis and those derived from the ORL analysis. Moreover the gas den- 
sity was gradually scaled up until the observed He 1 5876 A line flux 
was matched by the model. This approach, however did not succeed 
in reproducing the observed emission line spectrum, since the heat- 
ing, which, in the virtual absence of hydrogen, mainly comes from 
helium photoelectrons, was not adequate. This resulted in the knot 
being too cold (~300 K) to collisionally excite the optical and ul- 
traviolet heavy element forbidden lines observed in the spe ctrum. 
A sim ilar problem was encountered by lHarrington & Feibelmanl 
1 1984) when they tried to construct photoionization models of the 
hydrogen-deficient knots in Abell 30. This led them to hypothesize 
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Table 1. Model elemental abundances in the two assumed regions (dense, 
cold core and thin hotter outer region) of the knot J3 in Abell 30. Also listed, 
in columns 3 and 6, respectively, are empirical elemental abundances from 
the whole knot, derived from ORLs and from CELs Wesson et al. 1 2003). 
The abundances are given by number with respect to helium. 





Core 


ORLs 


Env. 




CELs 




H/He 


0.0250 


0.0850 


0.025 




0.085* 




C/He 


0.0200 


0.0389 


1.25-10" 


-3 


1.41-10" 


4 


N/He 


0.0222 


0.0229 


1.87-10" 


-3 


6.76-10" 


■5 


O/He 


0.0233 


0.1071 


3.50-10" 


-3 


1.78-10" 


■4 


Ne/He 


0.0200 


0.0831 


4.00-10" 


-3 


5.13-10" 


4 



* No es timate for the helium abundance was available to Wes son et alj 
1 2003) from CELs. H/He was assumed to be equal to the value they 
obtained from the ORL analysis. 



Table 2. Adopted parameters for the final Abell 30 knot J3 photoionization 
model 



log(L/L Q ) 


3.70 


R-core 


0.45-R fcnot 


T» 


130,000 K 


N He (core) 


10,400 cm" 3 


l°gg 


6.0 


N ffe (env) 


1680 cm" 3 




1.68-10 17 cm 


Pdlpg (core/env) 


0.077/0.107 


Rfcnot 


4.50-10 15 cm 


3d 


1.0-10 -6 cm 



Dknot- distance of the centre of the knot from the ionizing source 

Rfcnot : radius of the hydrogen-deficient knot 

R core : the radius of the thick, cold inner core of the knot 

N^ e (core): helium density by number [cm" 3 ] in the inner core of the knot 

Nue(env): helium density by number [cm" 3 ] in the outer region of the 

knot 

Pdlpg '■ dust-to-gas ratio by mass in the core and in the envelope of the knot 
a^ : radius of the dust grains 



a chemical inhomogeneity in the knot, as mentioned in SectionQ 
with a central region where the ratio of C, N and O relative to He is 
very high (~10%). However, as it was difficult to construct a one- 
dimensional photoionization model capable of self-consistently re- 
producing the postulated chemical inhomogeneity, Harrington & 
Feibelman could not test their hypothesis. 

With their result in mind, we constructed a second model, con- 
sisting of a spherical blob with two density/composition phases, 
where the core is assumed to be much denser and more metal- 
rich then the surrou nding envelope. The abundances derived by 
IWesson etafl f2003) from the CELs were taken as an initial guess 
for the model abundances in the envelope of the knot and the abun- 
dances from the ORLs were taken as a starting point for the abun- 
dances in the core. However, the He/H abundance ratio cannot be 
derived from the observed CELs. This was set to a starting value of 
0.5 by number, i.e. enhanced by a factor of 5 with respect to typical 
values for other PNe, the same factor by wh ich the C/H abundance 
appeared to be enhanced from lWesson et all l2003l) 's CEL analy- 
sis. The abundances were iteratively adjusted as we tried to fit the 
observed emission line spectrum. Ultimately, this two-phase model 
also failed to produce enough heating for the envelope of the knot, 
which, according to Wesson et al.'s analysis, has an electron tem- 
perature of Te(Om) = 16680 K. 



2.2 Dust heated models 

The failure of the first two models shows that photoionization of 
He and heavy elements alone cannot, in the virtual absence of 
hydrogen, provide enough heating to account for the collision- 
ally excited emission lines observed, so an extra source of energy 
must be responsible for the heating necessary to produce the ob- 
served CELs. Harrington & Feibelman jl984l) proposed that the ex- 
tra heating source might be provid ed by electron con duction from 
the shocked stellar wind. However, Harrington 1 1996) later argued 
that, while this is a possibility, the presence of even a weak mag- 
netic field could prevent such conduction and that a better candi- 
date for the additional energy source, needed to produce the ex- 
tra heating, is photoelectric emission from dust grains. This was 
first proposed by iBorkowski & Harrington! Il99lh in their analy- 
sis of the hydrogen-deficient planetary nebula IRAS 18333-2357. 
Evidenc e for the presence of dus t in Abell 30 is provided by the 
fact that lCohen&Barlowl ll974l) detected 10- and 20-/^m emis- 
sion from warm dust located w ithin a radius of 5.5 arcsec of the 
central star, while Icohen et alJ fl977) detected 1.6-3.5-/xm emis- 
sion from hot dust grains located in the same region. In addi- 
tion, the central star appears to be reddened, and, as discussed by 
lHarrington & Feibelmanl ll984l) . given the high galactic latitude 
of thi s object, the redd ening is likely to be intrinsic. A study by 
IBorkowski et alJ Il994l) . based on optical and near infrared images, 
confirmed the presence of dust in an expanding equatorial ring of 
H-poor gas around Abell 30; they detected no emission from the 
H-poor polar knots, concluding that there must be a dust deficiency 
there relative to the equatori al ring. Th is was confirmed by a more 
recent study bv lHarrington et ai]ll997l) . which showed the strength 
of the CELs relative to the He recombination lines to be much 
higher in the equatorial ring knots. They speculated that heating by 
photoionization might suffice in the polar knots, but that grain heat- 
ing had to be dominant in the equatorial knots. The results of this 
paper show that dust heating is in fact needed in the polar knots, in 
order to reproduce the correct line intensities for CELs. The iplica- 
tion of the dust deficiency of the polar knots to the equatorial knots 
and the results of the current work is that the dust-togas ratio in 
the equatorial knots must be much higher than that fou nd here for 
the polar knots. In the recent analysis of lWesson et al. 1 2003J), on 
the other hand, it was found that the extinction coefficient, c(Hf3), 
derived from the spectrum of the central star (c(H(3) = 0.60) and 
that derived for the knot J3 (c(H/3) = 0.55) are in good agreement. 
It is hard, however, to estimate the amount of dust actually locked 
into these knots, and, therefore, in this work, we treated the dust- 
to-gas ratio as a free parameter. A radiative transfer treatment for 
dust grains has not yet been implemented in the radiative transfer 
of MOCASSIN, although work is already in progress in this direc- 
tion. Nevertheless, in this work we estimated the heating provided 
by photoelectric emission from t he grains by using a strategy sim- 
ilar to that already used by Borkowski & Harrington 1 1991) for the 
modelling of the grain-heated, dusty planetary nebula in M22. An 
outline of the method used in this work is given in the next section. 

The final model consisted of a spherical knot composed of 
two phases: a dense inner nucleus (A^HcCcore) « 10,500 cm -3 ), 
with a radius equal to 0.45 of the total radius of the knot, -Rknot 
(where i?knot =4.50-10 lo cm), surrounded by a less dense outer re- 
gion (A^HcCenv) « 1700 cm -3 ), occupying the remaining volume 
of the knot. Although the central core region only occupies about 
9.1% of the total volume of the knot, its higher density and higher 
metal abundance means that a significant fraction of the total mass 
of the knot (43.9%) resides in this volume. As for the preliminary 
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modelling the starting abundances were taken to be the empirical 
values obtained from CELs and from ORLs for the envelope and 
the core regions of the knot, respectively. The abundnaces were 
then adjusted, together with the other model parameters, in an it- 
erative fashion in order to fit the the onbserved emission line spec- 
trum. The final elemental abundances adopted for the two regions 
are summarised in TableQ together with the empirical abundances 
derived by Wes son et al . (2003). The abundances for the outer re- 
gion and for the cold dense core, where most of the ORL emission 
is produced, are given by number with respect to helium (columns 4 
and 2, respectively, in Table0, since this is the most abundant el- 
ement. The abundanc es in the core are co mpared to the empirical 
ORL abundances of Wess on et alj 120031) . also given by number 
with respect to helium in columns 2 and 3 of Table Q Most of 
the CEL emission is believed to be produced in the hotter enve- 
lope of the knot and, therefore, the model abundances in this region 
(column 4 of Table0 are compared with the empirical CEL abun- 
dances (column 5 of Table0. It is clear from TableQ however, that 
the final model abundances are quite different from the starting val- 
ues which were taken to be equal to the values derived empirically. 
It should be noted at this point that, in Wesson et al.'s analysis, the 
abundances were derived without taking into account the possibil- 
ity of chemical and density inhomogeneities within the knot, which 
would give rise to regions with different ionization structures. This 
implies a certain degree of freedom in the choice of abundances for 
our photoionization model. In the light of this, the discrepancies 
existing between model and empirical abundance values are easily 
explained. We discuss in Section 4 the implications of our model 
abundances in the light of Wesson et al.'s empirical analysis. 

Table [2] lists the parameters used for the final simulation. 
We found that the ionization structure of this nebula was best re- 
produced by using a H-deficient model atmosphere for the ion- 
izing source, with a stellar luminosity of L* = 4965 Lq, an ef- 
fectice temperature of T, = 130,000 K and gravity defined by 
log g = 6.0. These parameters are in fair agreement with the val- 
ue s of T ef f = 110,000 K, log g = 6.0 and L„ = 5000 L derived 
by Ha rrington & Feibelman 1 19841) using the Zanstra method, and 
also with the results of Leuenhagen et al. (1993), who derived val- 
ues of T e fl = 1 15,000 K and log g = 5.5 from their spectral analysis 
of the central star of Abell 30. The hydrogen-deficient non-LTE 
stellar atmospheres used in this work were calculated by using the 
Tubingen NLTE Model Atmosphere Package, as discussed in Sec- 
tion l2.4l Both the dust-to-gas ratio and the grain size have a direct 
influence on the amount of heating produced by the photoelectric 
emission process. It is worth noting at this point that pd/ pg was 
treated as a free parameter, and different ratios were tried for the 
two regions of the knot, until the heating rate was high enough to 
reproduce the observed spectrum. 



2.3 Grain Photoemission 

Photoelectric emission from dust grains occurs when electrons are 
ejected from a grain following the absorption of a photon. The ra- 
tio between the number of photoelectrons emitted and the number 
of absorbed photons defines the photoelectric yield, Y . This, of 
course, depends on the grain species and size, since a photoelec- 
tron could be excited by a photon several layers from the surface 
on the grain. In small grains the probability of the photoelectron 
being excited near the surface is higher than in larger grains, and, 
therefore, there is also a higher probability of the photoelectron es- 
caping instead of being reabsorbed somewhere else inside the same 
grain. Experimental data exist for several materials (e.g. the data 



for graphite and vitreous carbon by iFeuerbacher & FittorJI 1 97 
but for the present calculations the following approximation for the 
photoelectric yield is used <Drainell978T) 



y(M = *.(i-£ 

Y(hv) = 



hv^B 
hv < B 



(1) 



where B is a parameter which is chosen to be some energy above 
the photoelectric emission threshold energy, or work function, 
E w , below which the yield becomes negligible, and Y^, is an- 
other parameter which depends on the material used and is de- 
termin ed experimentally. In the standard model define d by Draine 
( 1978 ), based on the res ults of (Feuerbacher & Fitton ( 1973) an d 
Feuerbacher et al. 1 1972), which is appropriate for organic, metal- 
lic and silicate grains, Yoo = 0.5 and B = 8 eV, where Draine's small 
grain size yield enhancement factors have been incorporated. 

Another important quantity which must be considered in es- 
timating the grain heating is the energy distribution of photo- 
electrons ejected from the dust grains. In this calculation the 
photoelectron distribution function, f(hu,E), where E is the 
photoelectron energy, is approx imated using the formulation of 
lBorkowskT& Harrinatoi] ll99lh for carbon grains, valid for pho- 
ton energies above 1 1 eV, 

E : E<2.7eV 

f(hv,E)cc [ E _ h7 + a s {E _ z7) s^ . £>2 . 7eV (2) 

where a = 0.059 (26.4 - hv) for hv < 26 .4 eV and a = for higher 
photon energies. For a particular photon energy, hv, this photo- 
electron energy distribution is truncat ed at an electron ener gy of 
hv — E w , where E w = 5eV iFeuerbach er & Fittonlll972Tl . For 
photon energies below 11 eV, equation [2] is scaled by replacing E 
by (11 — E w ) E I (hv — E w ). This scaling is necessary be- 
cause, from the experimental data, it is apparent that there is a shift 
in the energy at which the photoelectron energy distribution func- 
tion, /, peaks, from ~2.7 eV at high incident photon energies, to 
lower electron energies for hv < 1 1 eV. In practice, f(hv, E) is the 
fraction of all electrons ejected by a photon of frequency v which 
emerge with kinetic energies between E and E + dE. 

In this work only a qualitative description of the role of dust 
in the energy balance of the nebula is adopted. We use a simpli- 
fied approach, where grain charge and other factors are not taken 
into account, and the frequency-dependent contribution of the dust 
grains to the heating of the nebula, Gd(v), is estimated at each grid 
cell according to 



G d (v) 



iid 



(tto!) ■ F{v) ■ Y{v)- <E V > 



(3) 



where rid is the local dust density by number, a d is the radius of 
the dust grains assumed to be constant throughout the ionized re- 
gion, and F(v) is the photon flux. The quantity < E v > is the 
mean photoelectron energy, weighted by the photoelectron energy 
distribution function, f(hv, E), namely, 



< E v > = 



Jo 



f(hv,E)EdE 



r 

Jo 



f(hv, E) dE 



(4) 



The total dust contribution to the heating is then obtained by 
integration of equation [5] over all photon energies above the work 
function, E w = hvth, 



G d 



Gd(v) dv 



(5) 
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where Umax is the higher limit of the frequency range considered; 
in the MOCASSIN models of Abell 30, v ma x was set to 24 Ryd 
(« 326.4 eV). 

Finally, the total dust contribution to the heating, calculated at 
each grid cell, is added to the local thermal balance equations as 
an extra heating channel, so that the equation of thermal balance 
simply becomes 

Gphot + Gd — L c + L r + Lff (6) 

where G v hot is the heating due to the photoionization of the gas 
(mainly helium, in this case, but all other elements are also taken 
into account), L c is the contribution due to collisionally excited 
radiative cooling and L r and Lff are the contributions due recom- 
bination and free-free radiation, respectively. 



2.4 Synthetic ionizing spectrum of the central star of Abell 30 

For the central star it is more realistic to use NLTE model atmo- 
sphere fluxes instead of blackbody "spectra" or simple H+He com- 
posed models iRauch 2002). Consequently, we use synthetic spec- 
tra which have been calculated from line-blanketed NLTE model 
atmospheres (Werner 1986; Rauch 1997). These are plane-parallel 
and in hydrostatic and radiative equilibrium. 

iLeuenhaeen et alj ll993l) presented a spectral analysis of the 
central star of Abell 30. They found T cff = 115, 000 K, log g = 5.5, 
He:C:N:O=41:40:4:15 by mass, and log M/Mgyr -1 = -7.3. 
Wind effects may have an influence on the flux level in the UV 
lGableretalll989h . However, since the mass-loss rate is relatively 
low, these will be neglected in our analysis. 

We have calculated He+C+N+O models with 
T c fl — 110 — 150, 000 K, log g — 6 (cgs), and an abundance 
ratio of He:C:N:O=41:40:4:15 by mass. The statistics of the 
model atoms are summarized in Table|3| The observed nebular 
emission-line spectrum is best reproduced by our MOCASSIN 
photoionization models when a synthetic spectrum calculated 
from our stellar He+C+N+O model atmosphere grid, with 
T c g = 130, 000 K and logy = 6 is employed as the ionizing 
source. In Figure \\\ the fluxes obtained using a hydrogen-deficient 
model are compared to those obtained using a black body and a 
hydrogen-rich model with T e ff = 130, 000 K and logg = 6. 
A comparison between the nebular results obtained using a 
hydrogen-deficient and a hydrogen-rich model with normal cosmic 
abundances, for the same central star and nebular parameters is 
presented in Section l3~4l 



2.5 Modifications to the mocassin code 

Some modifications to the original version of the MOCASSIN code 
had to be made in order to include the dust photoelectric heating 
treatment described in Section l2~3l They mainly consisted of the 
addition of a dust heating term (Equation|5J to the thermal balance 
equation (Equation These modifications, however, will not be 
included in the public release version of MOCASSIN, which will 
consist of the pure photoionization version of the code only, since 
the intention is to enhance the code in the near future by including 
dust grains in the radiative transfer treatment. The analysis carried 
out in the current paper can, therefore, be considered to be a qual- 
itative investigation of the effects of dust photoelectric heating in 
hydrogen-deficient environments. 



Table 3. Statistics of model atoms used in the calculation of the NLTE stel- 
lar models. NLTE is the number of levels treated in NLTE, RBB is the 
number of line transitions. 194 additional levels are treated in LTE. 



ion NLTE RBB 



He I 


1 





Hell 


14 


78 


He III 


1 




C III 


3 


1 


C IV 


54 


288 


C V 


1 





N IV 


3 


1 


N V 


44 


150 


N VI 


1 





O V 


6 


3 


O VI 


43 


166 


O VII 


1 






Table 4. Mean temperatures (K) weighted by ionic species for knot J3 of 
Abell 30. For each element the upper row is for the optically thick inner 
core and the lower row is for the optically thin envelope of the knot. 



Ion 



El. 


I 


II 


III 


IV 


v 


VI 


H 


890 
17,270 


970 
17,320 










He 


500 
10,540 


950 
15,550 


1,430 
17,660 








C 


630 
11,600 


750 
11,820 


970 
15,830 


1,200 
17,500 


1,520 
17,780 




N 


550 
11,340 


630 
11,760 


960 
16,320 


1,200 
17,580 


1,510 
17,780 


1,520 
17,690 


O 


550 
11,310 


670 
11,760 


1,010 
16,460 


1,500 
17,690 


1,550 
17,800 


1,650 
17,690 


Ne 


500 
11,510 


690 
11,910 


1,040 
16,840 


1,510 
17,730 


1,550 
17,670 





3 ABELL 30 MODEL RESULTS 

3.1 The temperature and the ionization structure 

Table|4|lists mean temperatures (given to an accuracy of four sig- 
nificant figures), weighted by the ionic abundances, for the two 
regions of the knot. The upper entries in each row of Table |4| 
are for the inner core and the lower entries are for the surround- 
ing envelope. The value of T e (0 III) = 16,460 K obtained by the 
MOCASSIN model for the outer region of the kno t is in agree- 
ment w ith the empirical value derived from CELs bv lWesson et alJ 
(2003) of T e (0 III) = 16,680 K, as is also confirmed by the 
[O III] (4959 A + 5007 A)/4363 A model and empirical ratios that 
are in good agreement. The temperature corresponding to the pre- 
dicted [O III] (4959 A + 5007 A)/4363 A ratio of 52.4 is 17,000 K, 
which is about 500 K higher than the mean temperature weighted 
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Figure 1. Comparison of two synthetic NLTE model atmosphere fluxes calculated with T e tt = 130, 000 K, logg = 6, and different chemical com- 
position (Thin black line: H-rich model with normal cosmic abundances -labelled 'H-Ni'— ; grey line: He+C+N+O model with an abundance ratio of 
He:C:N:O=41:40:4:15 by mass). Prominent absorption thresholds in the He+C+N+O spectrum are marked. The dashed line (labelled 'BB') denotes the 
flux of a black body with T = 130.000K. 



Table 5. Fractional ionic abundances for knot J3 of Abell 30. For each el- 
ement the upper row is for the optically thick inner core and the lower row 
is for the optically thin envelope of the knot. Exponent powers of ten are 
given by the numbers in brackets. 



Ion 



El. 


I 


II 


III 


IV 


V 


VI 


H 


■ 19C-1) 
.77(-3) 


.981 
.999 










He 


.56(-l) 
.23(-3) 


.904 
.161 


■39(-l) 
.837 








C 


.66(-3) 
■27(-5) 


.112 
■71(-2) 


.800 
.139 


■87(-l) 
.530 


■24(-3) 
.324 




N 


.67(-3) 
■27(-5) 


.140 
.68(-2) 


.620 
.191 


.238 
.694 


.61(-3) 
.106 


.16(-6) 
■18(-2) 


O 


■ 82(-3) 
.16(-5) 


.158 
.66(-2) 


.820 
.278 


■21(-1) 
.627 


.42(-4) 
.89(-l) 


■19(-8) 
■17(-3) 


Ne 


.30(-2) 
.78(-5) 


.183 
■ 15(-2) 


.811 
.360 


.28(-2) 
.583 


.50(-5) 
■43(-l) 





by the O III abundance, as reported in Table [4] This can be ex- 
plained by the fact that there could be fluctuations in temperature in 
the O + region, in which case the empirical temperature derived by 
the [O III] diagnostic ratios would be weighted towards the higher 
limit of the range. 

The volume-averaged fractional ionic abundances for the two 
regions of the knot are listed in Table [5] where, once again, the 
upper entries for each element are for the inner core and the lower 
entries are for the surrounding envelope of the knot. 

It can be seen from Tables|4|and|5|that the model predicts very 
different temperature and ionization structures for the two regions 
of the knot. This has important implications for empirical deter- 
minations of abundances from spectroscopic observations, which 
generally assume constant ionization fractions throughout a region. 
In this light it is possible to understand the differences between the 
elemental abundances inferred using MOCASSIN models and those 
obtained by Wesson et al. (2003) using empirical methods. 

Finally, Figure |2| shows a slice plot of the temperature dis- 
tribution in the knot. In this figure the knot is oriented such that 
the stellar ionizing radiation is coming from the left-hand side. The 
large temperature variations existing in this model are illustrated by 
this plot, as well as the cooler shadow region created by the thick 
inner core. 
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Table 6. Dereddened observations and predicted emission lines fluxes for 
knot J3 of Abell 30. The line fluxes are given relative to He I 5876 A, on a 
scale where He I 5876 A = 1 . The flux of the He I 5876 A line is also given 



absolute units [10 ^^erg 


— 2 — 1 
cm s 


J- 






Line 




Predicted 




Observed 11 


Core 


Envel 


Total 




He I 5876/1Q- 15 er S 

cm - s 


7.07 


0.158 


7.23 


7.37 


He t 5876 


0.979 


0.021 


I 


I 


Hp tt Af.Sf, 

L1C 11 tUOU 


1 24 


869 


2 11 


2.50 


H/3 


062 


0023 


064 


080 


C II 4267 


347 


0002 


347 


0.452 


C nil 1 908 


0. 


1.52 


1.52 


2.8 


C ttt 4648 


0.044 


002 


046 


030 


C III 2297 


1.53 


0.016 


1.53 


0.960 


Civ 1550 6 


0. 


15.42 


15.42 


9.2 


[N II] 6584 


0. 


0.141 


0.141 


0.829 


N IV] 1486 


0. 


4.16 


2.62 


2.56 


[O II] 3727 


0.030 


0.120 


0.150 


1.66 


ii 4303 


0.030 


0. 


0.030 


0.094 


II 4097 c 


0.037 


0. 


0.037 


0.174 


On 4089 


0.072 


0. 


0.072 


0.060 


ii 4075 


0.158 


0. 


0.158 


0.149 


[0 m] 5007 


0. 


10.7 


10.7 


14.2 


[0 in] 4363 


0. 


0.275 


0.275 


0.275 


[Ne III] 3869 


0. 


7.57 


7.57 


2.85 


[Ne iv] 2423 


0. 


18.22 


18.22 


15.2 


[Ne v] 3426 


0. 


1.18 


1.18 


1.20 



a IWessonet"ai]|2003h . 

b Can be attenuated by dust absorption. 

c Observed line probably blended with the secondary Bowen fluorescence 
line from N III. 

3.2 The nebular emission line spectrum 

The adopted parameters for the models are those listed in Ta- 
bles 1 and 2. The fluxes of some important emission lines pre- 
dicted by MOCASSIN for the final model of knot J3 in Abell 30 are 
listed in Table |6] together with the corresponding observed values 
IWesson et all2003l) . The fluxes are given relative to He I 5876 A, 
since helium is the most abundant element in this object. The con- 
tributions from each of the two assumed regions of the knot are 
given separately in the table and, as expected, it is clear that in our 
model the ORLs are mainly emitted in the cold inner core of the 
knot, while the CELs are emitted in the hotter envelope surround- 
ing it. The predicted line fluxes are in satisfactory agreement with 
the observed values, although some discrepancies still remain, as 
discussed in the following paragraph. A dust to gas ratio in the en- 
velope of 0.107 by mass was sufficient to produce a heating rate 
adequate for the production of the observed CELs. In the central 
core the dust-to-gas ratio was set to 0.077 by mass, although, as 
will be discussed later, photoelectric emission from dust grains is 
not the dominant heating mechanism in this region. 

The largest discrepancies are observed for the [O II] 3727 A 
and [N II] 6584 A lines. Bo th of these li nes, however, may have 
a recombination component lRubirtll986J) . This component is ac- 
counted for by MOCASSIN in the case of [O II] 3727 A, but not in 
the case of [N II] 6584 A. The statistical equilibrium calculation 
routines were modified in order to self-consistently calculate the 
contribution of recombination to the population of each of the 5 lev- 
els considered in the 2+ ion. As the changes to the statistical equi- 



librium routines were made globally, the contribution of recombi- 
nation to any of the collisionally excited lines considered in a simu- 
lation can be taken into account, given that the required recombina- 
tion coefficients are known and have been added to the atomic data 
files. In our mod el we used the recombination coefficients given by 
iLiuet all f2000 ). their calculations however only covered the range 
of electron temperatures from 5,000 K to 20,000 K, while the elec- 
tron temperature in the core of the knot (where most of the recombi- 
nation lines are expected to be emitted) are in the range from 500 K 
to 2000 K. Low temperatures calculations for these recombination 
coefficients are not yet available, although work is being carried out 
at UCL by P. J. Storey. In the current work we have fitted a power 
law to the available data, in order to extrapolate the coefficients to 
the low temperatures required. The inclusion of radiative recombi- 
nation contributions to the final flux of [O II] 3727 A resulted in an 
increase of 25% (from a value of 0.120 to a value of 0.150 with 
respect to He 1 5876 A); this increase, although significant, is not 
sufficient to remedy the discrepancy with the observed value of 
1.66 with respect to He 15876 A. Due to the high electron densi- 
ties in the core region of the knot, collisional de-excitation plays 
an important part in depopulating the 2p 3 2 D level of 2+ hence 
suppressing the 2 Ds 3- 4 S3 doublet transition at 3726,29 A. The 

2 ' 2 2 

same is believed to occur in the case of the [N II] 6584,48 A transi- 
tions which are also underestimated in our models. An intermediate 
transition phase between the core and the envelope might provide 
the solution to this discrepancy; in our models we used a clear-cut, 
bi-density, bi-chemistry model, while a gradual decay in density 
and metal abundances might have been a more realistic description 
- the lower densities in such a region would lead to stronger CEL 
emission. 

There are also discrepancies between MOCASSIN'S prediction 
of the fluxes for some of the optical recombination lines of O II 
and the observations. However, O II 4303 A, 4097 A and 4089 A 
all come from the 3d 4 F level and therefore their relative strengths 
are only determined by the branching ratios, which are well known 
(e.g., see iLiu et al J 19951) . The observed relative intensities of these 
lines do not conform to the theoretical branching ratios. As indi- 
cated in Table [6] the O II 4097 A line should not be taken into 
account, since this may be blended with one of the secondary 
Bowen fluorescence lines of N III. Nevertheless, the relative fluxes 
of O II 4303 A and O II 4089 A are also in disag reement with 
the theoretical branching ratios. iTsamis et alj l2003h also noticed 
a similar effect in the spectra of five H II regions, namely, M 17, 
NGC 3576, 30 Doradus, LMC N11B and SMC N66. They sug- 
gested that at the low electron densities characteristic of those neb- 
ulae, a breakdown of the statistical equilibrium may occur among 
the fine-structure levels of the parent 3 Po,i,2, ground term of re- 
combining 2+ , causing the observed behaviour. However, in the 
case of the high electron density (7V C ~12,500 cm -3 ) in the core 
region of the knot J3 of Abell 30, this is not expected to occur. 

Table Q lists the diagnostic ratios for electron density and 
temperature predic ted by MOCASSIN and those observed by 
IWessonetail i2003). Good agreement is obtained for most line ra- 
tios. A lower ratio (-1 1.3%) was obtained for O II 4075 A/4089 A, 
indicating that the MOCASSIN model predicts a slightly lower tem- 
perature (1,150 K) for the inner core of the knot than the value 
implied by the observed ratio (2,800 K). This could be adjusted by 
using a slightly higher value for the dust-to-gas ratio in the core 
region of the knot. 
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Figure 2. Cross-section map of the electron temperature distribution in the 
MOCASSIN model of knot J3 of Abell 30. The slice plot is obtained on a 
plane through the centre of the knot and it is oriented such that the ionizing 
stellar field is incoming from the left-hand side. 



Table 7. Diagnostic ratios for electron temperature, T e , and electron den- 
sity, N e . Ratios predicted by MOCASSIN (column 3) and observed ratios 
(column 5). 



Ion 


Lines 


Ratio 


Ratio 




(A) 


(MOCASSIN) 


(Observed) 


N e 


[0 ii] 


3726/3729 


1.68 


1.73 


T e 


He I 


5876/4471 


2.80 


2.76 


He I 


6678/4471 


0.802 


0.777 


Oil 


4075/4089 


2.19 


2.47 


[O III] 


(5007+4959)/4363 


52.5 


54.5 



3.3 The thermal balance 

From the analysis of the heating and cooling contributions obtained 
at each location in our model grid, it appears that dust grains play 
the leading role in the heating of the knot in the envelope, dominat- 
ing over the heating produced by photoionization of helium, which 
is the most abundant element, and by other elements (including the 
small quantity of hydrogen remaining in the knot). This is very 
clearly visible in the left-hand panels of Figure|4] where the heat- 
ing and cooling contributions due to all the energy channels con- 
sidered in our simulation are plotted as a function of radius. These 
include cooling by collisionally excited emission lines, cooling by 
free-free and recombination continuum radiation and heating by 
photoionization and photoelectric emission from dust grains. The 
corresponding helium ionic fractions and the electron temperatures 
are plotted on the right-hand panels of Figure [4] also as a func- 
tion of radius. The plots correspond to three cuts through the knot, 
namely cut A (top panels), cut B (middle panels) and cut C (bottom 
panels). These are labelled in Figure|5| 

Heating from atomic photoionization is dominant in the inner 




Figure 3. Schematical representation of a cross-section through the centre 
of the knot. The three rays, labelled A, B and C, correspond to the posi- 
tions for which the model thermal contributions are plotted in Figure^ The 
stellar field is assumed to be incoming from the left-hand side. 



regions of the knot since the abundance of singly ionized helium 
(or neutral helium, in the case of cut A) is highest there, due to the 
sixfold density increase in this region. The shadowing effect of the 
dense core is clearly visible in cut A, causing heating by photoion- 
ization to decrease sharply already within the inner region of the 
knot. Heating by photoelectron emission from dust grain is higher 
in the envelope of the knot, where the radiation field is stronger. In 
fact, it is clear from Equation [3] that Gd(i^), the frequency depen- 
dent contribution to the heating due to dust grains, is proportional 
to the flux, F{y) and to the dust-to-gas ratio. As noted by the ref- 
eree, the drop of 30% in the value of the dust-to-gas ratio from the 
envelope to the core is not sufficient to explain the drop in the grain 
heating shown in Figure|4| so this must be mainly due to the drop 
in F(y), especially the attenuation of radiation at hu > 54.4 eV. 
In the shadowed region, F(v) is dominated by the softer diffuse 
component of the radiation field, the stellar photons having been 
mostly absorbed within the dense core. For this reason dust heating 
will make a smaller contribution in the shadowed region, as con- 
firmed by Figure|4| 

It is worth re-iterating at this point that the emission from the 
surrounding nebular material is not accounted for in our model. 
This would also contribute to the diffuse radiation field impinging 
onto the shadowed region. However, due to the very low surface 
brightess of the nebular material, its contribution to the difffuse ra- 
diation field in the knot region is expected to be very small. 



3.4 Dependence of the nebular models on the ionizing 

radiation spectrum 

As discussed in Section 2.2 we have used hydrogen-deficient syn- 
thetic spectra for the central star, calculated from line-blanketed 
non-LTE stellar atmospheres. The results of our nebular model 
(model A) are compared, in Table 8, to those obtained using a H- 
rich stellar atmosphere model with the same values of T e / / and 
log g (model B). In an attempt to isolate the effects that using a 
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Figure 4. Left panels: Heating and cooling contributions for knot J3 of Abell, plotted as a function of radius for three different cuts (from top to bottom: cuts A, 
B and C as shown in Figurel3l. (photo: photoelectric heating by atoms and ions; dust photoelectric heating by dust grains; coll: collisional cooling by ions;^: 
cooling by free-free radiation; rec: cooling by recombination). Right panels: Corresponding fractional ionic abundances for helium and electron temperatures, 
plotted as a function of radius. 
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Table 8. Comparison of results from photoionization models of knot J3 of 
Abell 30 with the dereddened observations using a non-LTE H-deficient 
model atmosphere (model A) and a non-LTE hydrogen-rich atmosphere 
(model B) as central ionizing star. The line intensities are given on a scale 
where I(He I 5876 A) = 1. The Observed column lists the dereddened ob- 
served intensities, while the A and B columns list the results obtained, re- 
spectively, using H-deficient (model A, as given already in Table 4) and 
H-rich (model B) atmospheres for the central ionizing source. 



Parameter 


Observed 


H-def 


H-rich 






star (A) 


star (B) 


Up i SS7ft/l 0~ l 5 /Wo/rm2/ c \ 

ne ijo/o/iu ^ei^/Liii /s ) 


7 37 


7 23 




He I 5876 


1. 


1. 


1. 


He II 4686 


2.50 


2.11 


2.12 


H I 4861 


0.080 


0.064 


0.066 


C II 4267 


0.452 


0.347 


0.261 


C hi] 1908 


2.8 


1.52 


0.958 


C in 4648 


0.030 


0.046 


0.025 


C ill 2297 


0.960 


1.53 


0.687 


Civ 1550 


9.2 


15.4 


8.02 


[N ii] 6584 


0.829 


0.141 


0.390 


N iv] 1486 


2.56 


2.62 


1.37 


[O II] 3727 


1.66 


0.150 


0.109 


II 4303 


0.094 


0.030 


0.022 


O II 4097 


0.174 


0.037 


0.028 


O ii 4089 


0.060 


0.072 


0.054 


O ii 4075 


0.149 


0.158 


0.103 


[O in] 5007 


14.2 


10.7 


5.32 


[O III] 4363 


0.349 


0.275 


0.061 


[Ne III] 3869 


2.85 


7.57 


1.87 


[Ne IV] 2423 


15.2 


18.22 


7.28 


[Ne v] 3426 


1.20 


1.18 


7.3 



different ionizing spectrum for the central source might have on 
the nebular emission line spectrum, we have kept all other nebu- 
lar parameters fixed (see Tables 1 and 2). The H-rich (H-Ni) atmo- 
sphere used for model B is one of a set of model calculated with line 
blanketed NLTE model atmospheres ( Wer neJl 1 986t lRauchlll99% . 
The nebular emission line fluxes are listed, with the dereddened 
observed fluxes being given in the second column and the model 
values in the following two columns. The fluxes are given relative 
to He I 5876 A (which is given in units of 10~ 15 ergs cm -2 s _1 in 
row 1 of Table|8j, on a scale where /(He I 5876 A) = 1. It is clear 
that the observed nebular spectrum is more closely reproduced by 
the photoionization model which used a hydrogen-deficient model 
atmosphere for the central star (model A). The crucial failure of 
model B is in the predicted [O III] 5007/4363 ratio of 87.2, versus 
model A's value of 38.9, which more closely reproduces the ob- 
served value of 40.7. The much higher ratio predicted by model B 
is indicative of the lower electron temperature obtained for the O III 
region. 



4 CONCLUSIONS 

We have constructed a photoionization model for the hydrogen- 
deficient knot J3 of Abell 30. This consisted of a spherical blob 
with a very dense metal rich core, surrounded by a layer of less 
dense gas, with somewhat less enhanced metal abundances. The 
envelope region composition might be formed by the mixing of 
the outer layers of the knot with the surrounding nebular gas. Our 



models found that chemical and density homogeneities alone were 
not sufficient to create enough heating to explain the collisionally 
excited lines (CELs) observed in the spectrum of the knot (see, for 
instance. lWesson et all2003l) . 

Following the hypothesis of Borkowski & Harrington 1 1991), 
heating by photoelectric emission from dust grains was introduced 
into the model. Although using a simplified approach, this work 
has shown that for a grain radius of 10 -6 cm a dust-to-gas ratio, 
Pd/pg, of 0.107 by mass is sufficient to provide enough heating to 
the envelope of the knot to fit the CEL spectrum. Although the final 
value used for our model is about 10 times larger than the value of 
0.01 -0.02, thought to be typical for the interstellar medium, this 
is not surprising, as the actual amount o f dust mixed with t he gas 
in thes e types of object is very uncertain. Borkowski & Harrington 
ll99lh inferred even larger dust-to-gas mass ratios, between 0.17 
and 0.62, from their study of the dusty H-deficient planetary nebula 
in M22. However, as discussed in Section 2.2, the dust-to-gas ratio 
in the equatorial k nots of Abell 30 is exp ected to be higher than 
in the polar knots ( Harring ton et alJ ll997). hence higher than the 
value of 0. 107 found here for the envelope the polar knot J3. 

Our final abundances derived for the two regions of the knot 
show core/envelope abundance ratios ranging from 5 for neon to 16 
for carbon; these are somewhat less dramatic than the ratios of 300- 
700 fo und from the ORL/CEL abundance analysis of Wess on et alJ 
(2003), but nevertheless follow the same trend. It should also be 
noticed, when comparing our abundances with those derived by 
IWesson et alJ 1200 3t) . that the empirical results they obtained were 
based on the assumption of a chemically/density homogeneous 
knot, with a constant ionization structure. However, this is not 
likely to be the case since it is clear from Table 7 that the ioniza- 
tion structure is very different in the inner and outer regions of the 
knot. In particular, the abundance of oxygen derived from empirical 
analysis is very susceptible to the empirical ionization correction 
used and would change considerably if a bi-phase empirical analy- 
sis were to be carried out. Nevertheless, in agreement with Wesson 
et al.'s empirical analysis, our models find both phases of the knot 
to be oxygen-rich, with C/O ratios of 0.86 and 0.36 in the core and 
the envelope of the knot, respectively. This is in contrast to theoret- 
ical predictions for PNe in the born-again scenario which predicts 
C/O ratios larger than unity llben & Renzini 1983; Herwig 2001). 
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